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In the Kuramoto model, a uniform distribution of the natural frequencies leads to a first-order 
(i.e., discontinuous) phase transition from incoherence to synchronization, at the critical coupling 
parameter K c . We obtain the asymptotic dependence of the order parameter above criticality: 
r — r c oc (K — K c ) 2 ^ 3 . For a finite population, we demonstrate that the population size N may be 
included into a self-consistency equation relating r and K in the synchronized state. We analyze the 
convergence to the thermodynamic limit of two alternative schemes to set the natural frequencies. 
Other frequency distributions different from the uniform one are also considered. 
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I. INTRODUCTION 

Synchronization is a universal phenomenon that plays 
an important role in all natural sciences as well as in tech- 
nology In particular, the synchronization of pop- 
ulations of globally coupled oscillators with distributed 
natural frequencies has been an object of study since very 
early times, mainly in a biological context 3]. Later, it 
has found application in other areas, such as Josephson 
junctions nanomechanics @, etc. When increasing 
the coupling parameter, these systems undergo transi- 
tions from a totally incoherent state to a partially co- 
herent state where part of the population becomes en- 
trained sharing the same frequency. Interestingly, there 
are several analogies with the phase transitions in statis- 
tical mechanics 6j. Thus, one may define an order pa- 
rameter, that usually grows continuously from zero (the 
incoherent state) when the coupling parameter exceeds 
a threshold value, analogously to a second-order phase 
transition. Nonetheless, in some situations 0)0,0) mu- 
tual cntrainment occurs in an abrupt way (a first-order 
phase transition) . After an infinitesimal variation of the 
coupling strength a macroscopic (i.e. order N) part of the 
population becomes synchronized. One may speculate 
that first-order phase transitions may be of interest for 
practical applications, if one pursues a system exhibiting 
an abrupt off/on switch. 

The most simple example of a first-order phase transi- 
tion is found in the Kuramoto model |lfj| when the natu- 
ral frequencies are uniformly distributed. In this case, it 
is known that all the population becomes synchronized 
m a sine le step We obtain here the asymptotic 

dependence of the order parameter after criticality, which 
exhibits a critical exponent 2/3. 

Still, one of the open problems in the Kuramoto model 
is to fully understand the finite- iV behavior. As we show 
below, the simplicity of the uniform frequency distribu- 
tion allows to cope with finite-size effects in an original 
way, providing analytic and numerical results. Some of 



these results apply to other frequency distributions with 
compact support. 

This paper is organized as follows. In Sec. [HI we 
present the Kuramoto model, and show some numeri- 
cal results that motivated this work. In Sec. IIIII an in- 
finite population is considered, finding the asymptotic 
dependence of the order parameter after criticality. In 
Sec. IIVI we study finite ensembles, finding: a) an N- 
dependent formula for the order parameter, b) different 
convergence rates to the thermodynamic limit for alter- 
native sampling schemes of the natural frequencies. Sec- 
tion is devoted to analyzing frequency distributions 
different from the uniform one. Finally in Sec. IVII the 
main conclusions of this work are summarized. 



II. THE KURAMOTO MODEL 

The Kuramoto model is probably the most studied 
model of synchronization in a population of oscillators 
with all-to-all coupling [l^. The state of each oscillator 
is described only by a phase variable (this stems from the 
the fact that, at small coupling, only the phase of a self- 
sustained oscillator is affected by the interaction). The 
phase 9j of each oscillator satisfies the following ODE: 

K N 

^ = ^- + ^E sin ^-^) (!) 

i=i 

where u>j are the natural frequencies, and K is the pa- 
rameter controlling the coupling strength. 

To quantify the state of synchronization, Kuramoto 
proposed to use a complex-valued quantity (so-called or- 
der parameter to emphasize the relation with phase tran- 
sitions): 

1 N 

r ^ = jjY,^ ei - ( 2 ) 

i=i 

It allows us to set the governing equation Q in the form: 
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Qj = ujj + Kr sin(V> — Oj). 



(3) 
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If the natural frequencies are distributed (i.e. ujj ^ 
u)ji), synchronization only appears above some coupling 
threshold. Here, we consider the case of evenly spaced 
natural frequencies: 
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-7+^(2j-l) 



(4) 



which, like all symmetric frequency distributions can be 
assumed centered at zero (by going into a rotating frame 
if necessary). Throughout this paper, the numerical inte- 
gration of the Kuramoto model [Eq. JTJ] is carried out by 
means of a fourth-order Runge-Kutta method with time 
step At = 0.1. 

Recently, Maistrenko et al. ^| studied the Kuramoto 
model with a small number of oscillators and natural fre- 
quencies distributed uniformly (or close to that). They 
found that the synchronized state robustly splits into 
several clusters with different average frequencies. It is 
shown there, and in Fig. ^ that for N = 3 and 5, the 
synchronized state splits directly into TV clusters. But, 
for values of N other than 3 and 5 the scenario is not so 
simple. As N increases, the number of splittings for go- 
ing from one to N clusters, increases as well. Hereafter, 
we denote by K s the coupling at the frequency-splitting 
from the synchronized state. In congruence with the ther- 
modynamic limit (see below) , all the splittings [including 
the first one at K — K S (N)\ must accumulate at K c as 
N — * oo. K c is the abrupt transition point for an infinite 
population. 



III. INFINITE POPULATION 

In this section, we briefly study the Kuramoto model 
for a uniform frequency distribution. Some of the formu- 
las will be later compared to those obtained for a finite 
population in Sec. II VI We also go one step further, and 
deduce an explicit formula, with critical exponent 2/3, 
for the dependence of the order parameter after critical- 
ity. 

In correspondence to the finite case, Eq. (J1J, we con- 
sider a uniform density of the natural frequencies: 



^ for M < 7 
for \oj\ > 7 



(5) 



We first note that due to the invariance under global 
rotation, for stationary solutions, we can set a vanishing 
phase for the order parameter in Eq. I|2|l. ib = 0. without 
lack of generality. Kuramoto's classical analysis gives the 
order parameter equation: 
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FIG. 1: Average frequencies (pj) as a function of the coupling 
for different population sizes. The natural frequencies are 
taken according to (0J with 7 = 1/2. By K 3 we denote the 
value of K for the first splitting bifurcation. The critical point 
K c in the thermodynamic limit (N — > 00) is located at the 
dotted line. 



In the totally locked regime, we obtain then: 
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(7) 
(8) 



cos8(ui) g(uj) dui . (6) 



Eq. © gives implicitcly the dependence of r on K. A 
solution exists only for Kr > 7. At the critical point 
K c r c — 7, the locked solution disappears with r c = ir/4. 
The corresponding value of the coupling is K c = 47/71-, 
that is precisely the value where the incoherent solution 
r = becomes unstable according to the classical re- 
sult 01 f° r a ll unimodal distributions: K c — 2/(ng(0)). 
Two remarks are in order. First, in contrast with strictly 
unimodal distributions, the transition is of first-order 
type: the order parameter "jumps" from zero to r c . Sec- 
ond, at K c all the population becomes entrained. This 
last remark is quite important because it simplifies both 
numerical and theoretical analyses. Also, note that when 
synchronized, the oscillators' phases are spanned along 
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FIG. 2: Log-log dependence of r on if in a neighborhood of 
criticality for 7 = 1/2. The circles correspond to numerical 
solutions of the self-consistency condition itS) . the solid line 
depicts Eq. I|17^ , and the dotted straight line arises taking the 
leading 5K 2 ^ 3 term only. 



an interval (of length tt at K = K c ). When K is in- 
creased r grows from r c to 1 in the K — > 00 limit. 

The first result of this paper is the dependence of r 
on K, just above the phase transition. First of all, we 
make a change of variables onto Eq. as usual (see 
e.g. 0,111): 



r = Kr 



cos 2 g(Kr sin#) dO. 



(9) 



where max (9 m i n ) is the phase of the oscillator with fre- 
quency 7 (—7). After an expansion above criticality: 

K = K C + SK (10) 
r = r c + 8r (11) 

flmaz = -e min =7T/2-Se. (12) 

and discarding the trivial incoherent solution r = 0, we 
get: 
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(13) 



An expansion of the sine function up to the cubic term 
yields: 
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(14) 
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Thus the problem reduces to finding 88, from Eq. ©: 

7 = (# c + (5#)(r c + <fr)sin(7r/2-£0) (15) 
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(16) 



Introducing this expression into Eq. I|14[l. a formula for 
the asymptotic dependence of <5r on is obtained: 



8r = 
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(17) 



The result is compared in Fig. [3] to the exact solution, 
arising from numerically solving Eq. (jSJ) . It confirms that 
the order parameter grows from r c with a power of K—K c 
with exponent 2/3. To our knowledge, this exponent and 
expression (|17|l have not been reported before. 



IV. FINITE POPULATION 

Finite-size effects in the Kuramoto model have been 
previously considered in the literature. Among the dif- 
ferent approaches, we may list: the investigation of the 
divergence of fluctuations around criticality [|| 0] , the 
observation of ephemeral coherent structures in the in- 
coherent state [l7j , and the reduction to a normal form 
in the case of identical natural frequencies with additive 
noise 0. Also, very recently, Mirollo and Strogatz 
have analyzed the (local) stability of the fully locked state 
for a finite population, finding that the locked solution 
is stable and disappears in a saddle-node bifurcation 
(as already observed in for small N). 

In this section, we show that for the uniform frequency 
distribution finite-size effects — on the order parameter 
and on the loss of the synchronization — can be studied 
in a novel way. 



A. Dependence of the order parameter on K 

For a finite population in the synchronized state (K > 
K s ), the order parameter is expressed [in correspondence 
with the integral form Q] by: 



1 N 

-Y 



1 - 



(18) 



where the order parameter r is a time independent quan- 
tity. We devote the following lines to deducing an N- 
dependent self-consistency equation [that reduces to (JHJ 
in the N — > 00 limit]. It is accurate provided that K is 
not too close to K s . 

According to the natural frequencies of the finite 
population are taken with step Aw = = 2^//N. 

Therefore, the order parameter H18|) is equivalent to a 
Riemann sum with constant step of the integral found in 
the continuum limit Q. The discrepancy between both 
cases is, at the leading order, dependent on /" (where 
f(u) = y/l - (cu/Kr) 2 ). Hence, for one Riemann box 
centered at lo^: 



f(uj)dtu = f(bJj)Auj + f ^ A^ 3 + 0{Au 5 ). (19) 



Therefore, for the finite N case we may approximate (jl8f) 
by: 



in u z , 
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In this expression, the sum may be approximated by its 
corresponding integral 
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f"{uj)duj + 0{N- 2 ) 
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And we obtain a self-consistency equation for a finite 
population of N oscillators: 



r= 1 -sh 
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With respect to the equation for the thermodynamic 
limit [Eq. ijHJl], there is an additional term depending 
on N. Plots of the numerical solutions for N = 10 
and 100 are shown in Fig. [3J Eq. (|22[) is deduced from 
the continuum equation, and therefore cannot intersect 
the line r — j/K, because several terms explode. Note 
that adding more terms of the series in Eqs. (|19I21(I does 
not overcome this problem. This suggests that, unfortu- 
nately, the behavior very close to K a cannot be deduced 
by simply manipulating the equations for an infinite pop- 
ulation. Nonetheless, out of that region, the solution of 
Q22JI reproduces the numerical results, even for such a rel- 
atively small number of oscillators as N = 10, Fig.^a). 
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B. Thermodynamic limit of the first 
frequency-splitting 

The arrangement of the natural frequencies in Eq. (|4"|) 
converges to the uniform frequency distribution (JSJ). But, 
as explained above, this limit is non-trivial: a first-order 
phase transition is substituted, when N becomes finite, 
by a set of frequency-splitting bifurcations accumulating 
at K c . Numerically, the study of these bifurcations is 
quite involved. Nonetheless, the point where, as K de- 
creases, the first splitting occurs {K = K s ), can be accu- 
rately computed. This is possible because above K s the 
system is in a fixed point state. Our simulations, Fig. 0] 
indicate that for the arrangement in (0J K s converges to 
K(T according to a power law: 

K c -K s (N)(xN-^. (23) 

Note that K s and K c are both proportional to 7, so /i ~ 
1.5 is independent of 7. 

The recipe followed to mimic the thermodynamic limit 
was to divide the frequency distribution g{uj) in N parts 
of equal area taking each Uj at the center of each block. 

Nonetheless, there are many (infinite in fact) possible 
discrete arrangements with the same continuum limit, 
but for arrangement fi is large enough to deduce the 



FIG. 3: Dependence of r on K for finite populations: (a) 
N = 10, and (b) N = 100; 7 = 1/2 in both cases. Values 
obtained from a direct computation of the Kuramoto model 
are marked by x . The first frequency splitting is observed 
at Ks (for K < K s , r is no longer constant in time). The 
solid line is obtained solving Eq. I22H numerically. It matches 
the computed values, but fails when approaching the line r = 
y/K. As a reference, the solution for an infinite population 
is shown with circles [after numerically solving Eq. P)]: the 
critical point (K c ,r c ) is marked with a • symbol. 

decay for other distributions, by just considering the vari- 
ation of the effective 7. For instance, if the natural fre- 
quencies are taken as in 0, [2(| : 

Wi = -7+^rj-2Cj-l) (24) 

one observes that K s converges to K c + with a power-law, 
see o's in Fig. 0] but more slowly than for the arrange- 
ment in Q. As the extrema of (|24|) are fixed, a change in 
N varies the effective width of the equivalent continuous 
distribution: 7 e // = 7 + 7/(iV— 1). Hence, from 123|) we 
get: 

KJN) = ^it -aN'" ^K c + -^ + 0(N-") (25) 

which agrees with the observed result (see the dashed line 
in Fig. EI). 
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FIG. 4: Log-log dependence of the distance from K s to K c 
on the population size (7 = 1/2). Squares and diamonds 
correspond to different arrangements of the natural frequen- 
cies, Eqs. Q and 1241 . respectively. In the first case, the last 
two decades were fitted with a straight line, yielding a slope 
—fM = —1.502. The dashed straight line log 10 (47/7r) — x arises 
from Eq. ((23. 
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TABLE I: Three frequency distributions considered in Sec.lVl 
V. OTHER FREQUENCY DISTRIBUTIONS 



For both schemes and the three frequency distribu- 
tions in Table CI the approach of the first frequency split- 
ting to the thermodynamic limit satisfies a power law: 
— K S (N)\ oc N~^, as occurred for the uniform dis- 
tribution (K%° — K c in this case). For sampling (ii) 
we find that the value of the exponent is always fj, rs 1. 
However for (i), different exponents arise, in contrast to 
/i w 1.5 obtained for the uniform distribution: \l « 0.5 
for triangular and parabolic distributions, and fj, ~ 1 for 
the hat-shaped one. We have checked that this exponent 
arises for other distributions with an abrupt boundary 
(<7(±7) > 0). 

Another interesting power law is the shift of the or- 
der parameter in the synchronized state: \r(K,N) — 
r(K,N — 00) I ~ N~ y . For sampling scheme (i) theo- 
retical results may be obtained, using again arguments 
based on the Riemann sum. If g (±7) > — e.g. uni- 
form [see Eq. 1)22(1 ] or hat-shaped distributions — one may 
obtain a formal solution (23 that yields v = 2. How- 
ever for distributions that approach zero at u = ±7, the 
"Riemann-sum approach" is not valid due to divergences 
at ui = ±7. One must, therefore, analyze these points 
separately. In particular, one obtains v = 3/2 for lin- 
early decaying distributions [g(w — > ±7) ~ (7 =p u/), e.g. 
parabolic and triangular distributions]. Also, our simula- 
tions indicate that for sampling (ii) irrespective of 
the frequency distribution (for the uniform distribution 
v = 1 is straightforward due to the 0(N~ r ) effective shift 
of 7). 

Finally, we note that a recipe similar to (i) consisting 
in taking the frequencies at the median (instead of the 
center) of each block [24| exhibits the same exponents 
jtx, v that scheme (i). 



In this section, we briefly discuss the extension of the 
previous results to other frequency distributions sup- 
ported on a finite interval [—7,7]- We first note that for 
non-uniform distributions the lost of complete synchro- 
nization and the critical point where the incoherent solu- 
tion becomes unstable do not coincide: K%° = K S (N 
00) > K c . In other words, there exists an intermediate 
range of partial entrainment in which one part of the pop- 
ulation is synchronized whereas the remaining oscillators 
drift. 

As model distributions we considered three unimodal 
distributions listed in Table |U As for the uniform dis- 
tribution, the desynchronization point may be com- 
puted analytically using Eq. 0. 

We focused on two simple sampling schemes to dis- 
crctizc distributions supported on a bounded interval: 

(i) g{w)du = 2 1 /N 1 Q l = = £d£*±i. 

(ii) J^ +1 g{u)dw = 2j/(N - l),w N = -m = 7. 

Applying (i) and (ii) to a uniform distribution one gets 
the arrangements in Eqs. J3J and 1)24(1 . respectively. 



VI. CONCLUSIONS 

In the present paper, the first-order phase transition 
arising when imposing a uniform frequency distribution 
on the Kuramoto model has been studied. In the case of 
an infinite population, we have found an explicit asymp- 
totic dependence of the order parameter after criticality, 
Eq. (HI). 

For a finite population, our first conclusion is that 
in contrast to strictly unimodal distributions (Cauchy, 
Gaussian, parabolic,...) of the natural frequencies, which 
exhibit transitions of second-order type, the thermody- 
namic limit is non-trivial for a uniform distribution. In 
the finite- N case, the synchronized state does not split 
directly into N clusters, but through a cascade of fre- 
quency splittings. To be congruent with the first-order 
phase transition predicted in the thermodynamic limit 
all the splittings must accumulate at K c as N — > 00. 

The dependence of the order parameter r on the cou- 
pling K, has been expressed in an easy-to-compute for- 
mula, Eq. 1(22(1. In this formula, the population size N 
enters explicitly, and it allows, except very close to the 
desynchronization point, an accurate computation of the 
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order parameter, even for small N. 

Two sampling schemes to set the natural frequencies 
have been compared. The scheme we propose in Q con- 
verges to the thermodynamic limit faster than another 
used in the literature [Eq. I|24l) ]. The comparison is based 
on the different exponent of the power-law convergence 
for the point of the first frequency splitting (// ss 1.5 vs. 
H = 1) and the shift of the order parameter (v — 2 vs. 
v = 1). Other frequency distributions (with compact 
support) different from the uniform one have been dis- 
cussed. Among the infinite possible sampling schemes, 
those studied here appear the most natural ones to us. 
Nonetheless, further investigation is needed to assess the 
existence of an optimal sampling scheme to mimic the 
thermodynamic limit with a finite population. 

In spite of the lack of structural stability under per- 



turbations of the uniform frequency distribution (in the 
thermodynamic limit, but not in the finite case as proved 
in 0]), the results in this paper could be useful in order 
to understand more complicated schemes, like oscillator 
networks |2l|. The use of a frequency distribution with 
an abrupt transition seems more suited to better resolve 
critical points in this kind of systems. 
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